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Qs, ■ ABSTRACT 



Aims. In the present work we study the behavior of the magnetic field and the line activity indicators of the single 
late-type giant j3 Ceti. Using spectropolarimetric data, we aim to reconstruct the magnetic field structure on the star's 
■ surface and to present the first magnetic maps for j3 Ceti. 

Methods. The data were obtained using two spectropolarimeters - Narval at the Bernard Lyot Telescope, Pic du Midi, 
' France, and ESPaDOnS at CFHT, Hawaii. Thirty-eight circularly-polarized spectra have been collected in the period 

June 2010 - January 2012. The Least Square Deconvolution method was applied for extracting high signal-to-noise 
JL ratio line profiles, from which we measure the surface- averaged longitudinal magnetic field B; . Chromospheric activity 

indicators Call K, Ha, Call IR (854.2 nm) and radial velocity were simultaneously measured and their variability was 
analysed together with the behavior of B; . The Zeeman Doppler Imaging (ZDI) inversion technique was employed for 
reconstruction of the large-scale magnetic field and two magnetic maps of /3 Ceti are presented for two periods (June 
2010 - December 2010 and June 2011 - January 2012). 

Results. B; remains of positive polarity for the whole observational period and shows significant variations in the 
£SJ ' interval 0.1 - 8.2 G. The behavior of the line activity indicators is in good agreement with the B; variations. Searching 

, for periodic signals in the Stokes V time-series, we found a possible rotation period of 215 days. The two ZDI maps show 

£SJ ■ a mainly axisymmetric and poloidal magnetic topology and a simple surface magnetic field configuration dominated by 

^\ ' a dipole. Little evolution is observed between the two maps, in spite of a 1 yr interval between both subsets. We also 

' use state-of-the-art stellar evolution models to constrain the evolutionary status of /3 Ceti. We derive a mass of 3.5 Mq 

and propose that this star is already in the central-helium burning phase. 

Conclusions. Taking into account all our results and the evolutionary status of the star, we suggest that dynamo action 
alone may not be efficient enough to account for the high magnetic activity of f3 Ceti. As an alternate option, we propose 
that it may be an Ap star descendant presently undergoing central helium-burning and still exhibiting a remnant of 
the Ap star magnetic field. 
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?H ■ 1. Introduction main sequence, large-scale magnetic fields of cool evolved 
. y.' stars may be the product of a global dynamo, triggered 
Studying magnetic activity in cool single giants is crucial by the i nte rplay between internal differential rotation and 
for understanding the role of magnetic fields at the late convection in the stellar envelope (Konstantinova-Antova 
stages of stellar evolution. By analogy with cool stars on the ct al 2 010; 2012). Their magnetic fields may also be a 
fossil remnant inherited from previous evolutionary stages, 



Send offprint requests to: S. Tsvetkova still frozen in stable internal layers of the star and either 

* Based on observations obtained at the Bernard Lyot pervading convective layers up to the surface or interact- 

Telescope (TBL, Pic du Midi, France) of the Midi-Pyrenees hig with a global dynamo (Auriere et al. 2008; 2012). This 

Observatory which is operated by the Institut National des second scenario may be of particular importance for those 

Sciences de l'Univers of the Centre National de la Recherche cool giants which are the descendants of Ap/Bp stars for 

Scientifique of France and Universite de Toulouse, and at the which the fossil field hypothesis is generally preferred to 

Canada-France-Hawaii Telescope (CFHT) which is operated by acc0 unt for their strong and topological^ simple surface 
the National Research Council of Canada, the Institut National 
des Sciences de l'Univers of the Centre National de la Recherche 
Scientifique of France, and the University of Hawaii. 



magnetic fields. A third possible contributor to the surface 
magnetism of cool giants may be a highly-tangled mag- 
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netic pattern originating from a small-scale dynamo. While 
local dynamos are probably active in all stars with convec- 
tive zones, their direct observational signatures may only be 
detectable in the most evolved objects (cool supcrgiants), 
for which the larger pressure scale height results in much 
larger convective cells (Schwarzschild 1975). This last type 
of magnetism is presumably responsible for the weak mag- 
netic field detected at the surface of the slowly-rotating su- 
pergiant Betelgeuse (Auriere et al. 2010, Petit et al. 2013). 

Observations are obviously essential to decide between 
these different scenarios (that may act independently or 
together, depending on the object) or unveil other possi- 
ble physical explanations for the manifestations of atmo- 
spheric activity recorded in cool, evolved stars. Recent re- 
views of different types of active stars, including cool giants, 
have been presented by Berdyugina (2005) and Donati & 
Landstreet (2009). In addition to well-known and some- 
what atypical evolved objects (e.g. components of RS CVn 
systems and FK Com stars), weak surface magnetic fields, 
ranging from a few Gauss to a few tens of gauss, are now de- 
tected in more common, single G-K giants (Konstantinova- 
Antova et al. 2008; 2009; Auriere et al. 2009). Similar de- 
tections were also reported for a few M giants, e.g. EK Boo 
(of M5 spectral type, Konstantinova-Antova et al. 2010). 
The surface magnetic topology, and sometimes its tempo- 
ral evolution, has been reconstructed for a small number 
of cool giants - the Ap star descendant EK Eri (Auriere 
et al. 2011); the fast rotating effectively single giant V390 
Aur (Konstantinova-Antova et al. 2012); the early K-type 
primary giant of the RS CVn binary IM Peg (Berdyugina 
& Marsden 2006); and the FK Com star HD 199178 (Petit 
et al. 2004). 

P Ceti (HD 4128, HR 188, HIP 3419) is one of the mag- 
netic single late-type giants from the sample of Auriere et 
al. (2009). It is a star of spectral class K0 III with V = 2.04 
mag and B - V = 1.02 mag. Its value of V - R = 0.72 
mag places it to the left of the coronal dividing line in the 
Hertzsprung-Russcl diagram (hereafter HRD) where stars 
are supposed to have a chromosphere, a transition region 
and a corona. That dividing line was proposed by Linsky 
& Haisch (1979) near V — R = 0.80 mag and later con- 
firmed by the study of Simon et al. (1982) with a larger 
sample of stars. In the solar neighborhood (d < 30 pc), f3 
Ceti is the single giant star with the highest X-ray lumi- 
nosity \ogL x = 30.2 erg/s (Maggio et al. 1998; Hunsch et 
al. 1996), calculated assuming a distance of d = 29.5 pc 
based on the Hipparcos parallax (van Leeuwen 2007). This 
high X-ray emission reaches a level comparable to that of 
Capella (a member of a long period RS CVn association) 
and 9 1 Tau (K0 III, in the Hyades open cluster) and reveals 
an extended corona. An atmospheric model of /3 Ceti was 
performed by Eriksson et al. (1983), suggesting the exis- 
tence of coronal loops. Ayres et al. (2001) reported about a 
series of striking coronal flare events observed with EUVE 
during a period of 34 days starting on 1 Aug 2000, confirm- 
ing the strong activity of this single giant. 

To study the magnetic field and activity of j3 Ceti, we 
have collected high-resolution spectropolarimctric data in 
the period June 2010 - January 2012. We present here 
the first Zeeman Doppler Imaging (ZDI) maps of this star. 
In Section 2 we describe the observations and data reduc- 
tion. Section 3 contains the results of our measurements. In 
Section 4, we derive the mass and discuss the evolutionary 
status of f3 Ceti. We then propose in Section 5 an inter- 



pretation of the reconstructed large-scale surface magnetic 
topology of /3 Ceti and the variability of activity indicators. 
Finally, Section 6 contains our concluding remarks. 

2. Observations and data reduction 

Observational data were obtained with two twin fiber-fed 
echelle spectropolarimeters - Narval (Auriere 2003) which 
operates at the 2-m Bernard Lyot Telescope (TBL) at Pic 
du Midi Observatory, France, and ESPaDOnS (Donati et al. 
2006 a) which operates at the 3.6-m Canada-France-Hawaii 
Telescope (CFHT) of Mauna Kea Observatory, Hawaii. In 
polarimetric mode, both have a spectral resolution of about 
65 000 and a nearly-continuous spectrum coverage from the 
near-ultraviolet (at about 370 nm) to the near-infrared do- 
main (at 1050 nm) in a single exposure, with 40 orders 
aligned on the CCD frame by 2 cross-disperser prisms. 
Stokes I (unpolarised light) and Stokes V (circular polariza- 
tion) parameters are simultaneously obtained by four sub- 
exposures between which the retarders - Fresnel rhombs - 
are rotated in order to exchange the beams in the instru- 
ment and to reduce spurious polarization signatures (Semel 
et al. 1993). 

38 spectra have been collected for (3 Ceti in the pe- 
riod June 2010 - January 2012 (Table [TJ. First, reduced 
spectra were extracted using the automatic reduction soft- 
ware LibrcEsprit, developed for Narval and ESPaDOnS. A 
detailed description of the algorithm implemented in this 
software can be found in Donati et al. (1997). As a sec- 
ond step, the Least Squares Deconvolution (LSD) multi-line 
technique (Donati et al. 1997) was applied to all observa- 
tions. This widely-used cross-correlation technique enables 
the averaging of several thousand absorption atomic lines 
recorded throughout the echelle spectrum, generating a sin- 
gle Stokes I and V line profile. We employed here a line 
mask calculated for an effective temperature T e // = 5000 
K, log<7 = 3.0 and a microturbulence of 2 km/s, consistent 
with the literature data for j3 Ceti (Thevenin 1998; Hekker 
& Melendez 2007; Massarotti et al. 2008) and containing 
about 12 700 spectral lines. By doing so, the signal-to-noisc 
ratio (S/N) is increased to the point where weak polar- 
ized Zeeman signatures can be detected in all 38 spectra, 
thanks to a final S/N of up to 56 000. Actually, we have two 
more spectra of j3 Ceti taken in September 15, 2010, and 
November 20, 2010, but both were excluded from our anal- 
ysis due to low S/N values; measurements extracted from 
these low-quality spectra are not included in the tables and 
figures. 

As a first estimate of the magnetic field strength, we 
computed the longitudinal component of the magnetic field 
(Bj), using the first-order moment method (Rees & Semel 
1979; Donati et al. 1997; Wade et al. 2000). According to 
the multi-line model, we assumed here that the LSD line 
profiles possess a mean wavelength of 597 nm and a Lande 
factor of 1.28. The first moment was computed between ra- 
dial velocity boundaries set to encompass the whole velocity 
span of Stokes V signatures. 

We also took advantage of several classical activity prox- 
ies showing up in the spectra, by measuring the relative 
intensity relative to the continuum (R c ) for the line activ- 
ity indicators Ha at 656.3 nm and the Call IR triplet at 
854.2 nm. For the Call K line, we considered the inten- 
sity ratio between the line emission core and the intensity 
at 395 nm. We also measured the radial velocity (RV) of 
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Tabic 1: Measured longitudinal magnetic field, activity indicators, radial velocity and their uncertainties for /3 Ceti. In the 
first column capital letters N and E stand for Narval (TBL, Pic du Midi, France) and ESPaDOnS (CFHT), respectively. 



Inst. Date HJD Rot, phase Call K a Ha a Call IR a B,, [G] a, [G] RV [km.s' 1 ] 



E 


19 Jun 10 


2455368.123 


0.891 





.191 


0.005 





.211 


0.002 


0. 


134 


0.003 


4.4 


0.4 


13.36 


E 


21 Jun 10 


2455370.134 


0.901 


0. 


194 


0.004 


0. 


214 


0.001 


0. 


135 


0.003 


3.9 


0.4 


13.36 


E 


16 Jul 10 


2455395.140 


0.017 


0. 


210 


0.004 


0. 


220 


0.001 


0. 


141 


0.001 


7.4 


0.3 


13.36 


E 


17 Jul 10 


2455396.132 


0.022 





216 


0.005 





217 


0.001 


0. 


144 


0.002 


8.1 


0.4 


13.38 


E 


25 Jul 10 


2455404.138 


0.059 





.207 


0.012 





215 


0.001 


0. 


146 


0.000 


8.2 


0.9 


13.47 


E 


04 Aug 10 


2455414.006 


0.105 


0. 


.218 


0.005 


0. 


213 


0.001 


0. 


147 


0.001 


7.2 


0.4 


13.45 


N 


06 Aug 10 


2455415.633 


0.112 


0. 


197 


0.007 


0. 


211 


0.002 


0. 


149 


0.001 


7.0 


0.5 


13.44 


N 


16 Aug 10 


2455425.685 


0.159 





.201 


0.004 





.221 


0.002 


0. 


148 


0.003 


4.5 


0.5 


13.43 


N 


02 Sep 10 


2455442.699 


0.238 


0. 


189 


0.007 


0. 


218 


0.001 


0. 


147 


0.001 


2.8 


0.4 


13.36 


N 


19 Sep 10 


2455459.473 


0.316 





.187 


0.005 





.211 


0.002 


0. 


142 


0.001 


1.2 


0.4 


13.36 


N 


26 Sep 10 


2455466.500 


0.349 





198 


0.006 





215 


0.001 


0. 


147 


0.002 


0.1 


0.4 


13.35 


N 


05 Oct 10 


2455475.518 


0.391 





212 


0.007 





215 


0.002 


0. 


155 


0.002 


1.8 


0.7 


13.43 


N 


13 Oct 10 


2455483.475 


0.428 


0. 


230 


0.005 


0. 


213 


0.001 


0. 


165 


0.001 


3.6 


0.4 


13.44 


E 


15 Oct 10 


2455485.832 


0.439 





216 


0.003 





.223 


0.002 


0. 


155 


0.002 


3.9 


0.4 


13.43 


E 


17 Oct 10 


2455486.863 


0.444 


0. 


222 


0.006 


0. 


217 


0.001 


0. 


154 


0.001 


3.8 


0.3 


13.42 


N 


20 Oct 10 


2455490.471 


0.460 





211 


0.016 





213 


0.001 


0. 


158 


0.001 


5.1 


0.8 


13.46 


N 


12 Nov 10 


2455513.360 


0.567 





199 


0.004 





212 


0.002 


0. 


149 


0.004 


4.9 


0.7 


13.40 


E 


15 Nov 10 


2455516.831 


0.583 





.202 


0.003 





.222 


0.002 


0. 


142 


0.003 


4.0 


0.4 


13.35 


E 


20 Nov 10 


2455521.798 


0.606 


0. 


197 


0.007 


0. 


212 


0.002 


0. 


137 


0.001 


4.2 


0.5 


13.31 


E 


21 Nov 10 


2455522.718 


0.610 





.195 


0.006 


0. 


212 


0.004 


0. 


137 


0.004 


3.9 


0.4 


13.34 


N 


26 Nov 10 


2455527.368 


0.632 





.185 


0.004 





212 


0.001 


0. 


136 


0.003 


4.5 


0.4 


13.26 


E 


27 Nov 10 


2455528.855 


0.639 





196 


0.006 





212 


0.002 


0. 


137 


0.002 


3.9 


0.4 


13.32 


N 


04 Dec 10 


2455535.349 


0.669 





.175 


0.009 





207 


0.001 


0. 


138 


0.003 


4.8 


0.6 


13.35 


N 


12 Dec 10 


2455543.383 


0.706 





175 


0.010 





211 


0.001 


0. 


139 


0.002 


5.0 


0.8 


13.32 


N 


14 Dec 10 


2455545.292 


0.715 





167 


0.013 


0. 


208 


0.002 


0. 


147 


0.003 


5.3 


1.1 


13.24 


E 


18 Jun 11 


2455732.097 


0.584 





211 


0.016 





.204 


0.003 


0. 


136 


0.004 


4.3 


0.8 


13.41 


E 


21 Jun 11 


2455735.130 


0.598 


0. 


213 


0.007 





205 


0.002 


0. 


135 


0.003 


4.1 


0.4 


13.41 


E 


07 Jul 11 


2455751.132 


0.673 





,227 


0.010 


0. 


208 


0.001 


0. 


145 


0.004 


4.7 


0.6 


13.42 


E 


12 Jul 11 


2455756.134 


0.696 





.225 


0.009 





.208 


0.003 


0. 


145 


0.005 


4.2 


0.5 


13.35 


E 


14 Jul 11 


2455758.138 


0.705 





.226 


0.007 





214 


0.003 


0. 


145 


0.005 


4.6 


0.4 


13.39 


N 


24 Sep 11 


2455829.535 


0.037 





205 


0.005 


0. 


222 


0.002 


0. 


157 


0.002 


6.7 


0.5 


13.39 


N 


10 Oct 11 


2455845.407 


0.111 





206 


0.008 





223 


0.002 


0. 


151 


0.002 


6.0 


0.6 


13.38 


E 


30 Oct 11 


2455865.444 


0.204 





.183 


0.006 





.208 


0.005 


0. 


141 


0.002 


3.7 


0.5 


13.38 


N 


16 Nov 11 


2455882.425 


0.283 





.180 


0.009 





211 


0.004 


0. 


143 


0.005 


4.3 


0.4 


13.34 


N 


27 Nov 11 


2455893.309 


0.334 


0. 


191 


0.006 


0. 


216 


0.006 


0. 


147 


0.005 


2.7 


0.4 


13.35 


N 


08 Dec 11 


2455904.316 


0.385 





203 


0.006 





216 


0.002 


0. 


154 


0.005 


1.3 


0.4 


13.35 


N 


09 Jan 12 


2455936.252 


0.534 





209 


0.005 





216 


0.002 


0. 


159 


0.005 


4.0 


0.5 


13.25 


N 


22 Jan 12 


2455949.279 


0.594 





.184 


0.004 





.213 


0.001 


0. 


145 


0.003 


5.2 


0.4 


13.32 



the LSD Stokes I profiles, using the ccntroid of a simple 
Gaussian fit. The RV stability of Narval and ESPaDOnS 
was estimated by Moutou et al. (2007) to be of the order of 
30 m/s. The absolute accuracy of individual measurements 
relative to the local standard of rest is about 1 km/s with 
both spectropolarimetcrs (Auricrc et al. 2008). 

We exploited the rotational modulation of Stokes V sig- 
natures to reconstruct the magnetic topology of the star 
by means of the Zeeman Dopplcr Imaging tomographic 
method (ZDI; Semcl 1989; Donati & Brown 1997; Donati 
et al. 2006 b). We used an implementation of this algo- 
rithm where the surface vectorial magnetic field is projected 
onto a spherical harmonics frame, allowing us to easily dis- 
tinguish between the poloidal and toroidal components of 
the surface magnetic geometry (Donati et al. 2006 b). This 
method performs iterative adjustment of a simulated time- 
series of Stokes V profiles and the observed time-series. The 
synthetic Stokes profiles are calculated from an artificial 
star whose surface is divided into a grid of 2 000 rectan- 
gular pixels of roughly similar area. Each surface pixel is 
associated with a local Stokes I and V profile. In this simple 



model, the local synthetic Stokes I line profile is assumed to 
possess a Gaussian shape, with a depth and width adjusted 
to achieve the best fit between synthetic and observed line 
profiles. Assuming a given magnetic field strength and ori- 
entation for each pixel, local Stokes V profiles are calculated 
under the weak-field assumption (Morin et al. 2008; Petit 
et al. 2010; Morgenthalcr et al. 2012). The rotation period, 
inclination angle and v sin i adopted for our model will be 
discussed in Sect. 13.21 The linear limb darkening coefficient 
is set to 0.75, in agreement with Claret & Bloemen (2011). 
We limited the spherical harmonics expansion to I < 10 
since no improvement is found for the fits between modeled 
and observed LSD profiles for higher values of I. 

3. Results 

3.1. Longitudinal magnetic field, line activity indicators and 
radial velocity 

All Stokes V profiles of Ceti have a simple shape with 
a peak-to-peak amplitude varying from 9.7 x 10 _5 / c to 
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Fig. 1: Variations of B; with activity indicators Call K (upper left), Call IR (upper right), Ha (lower left) and radial 
velocity (lower right) in the period June 2010 - January 2012. 



6.4 x 10~ 4 / c in 2010 and from 1.3 x 10- 4 / c to 6 x 10~ 4 / c in 
2011/2012. All line profiles feature a positive blue lobe and 
negative red lobe, which corresponds to a field of positive 
polarity. Accordingly, B; estimates remain of positive sign 
for the whole observational period (see Table Q] and Fig. [I}. 
Thanks to a typical error bar as low as 0.5 G, significant 
variations are observed in the interval 0.1 - 8.2 G in 2010 
and 1.3 - 6.7 G in 2011/2012. The smaller amplitude of vari- 
ations observed in 2011/2012 may be plausibly explained 
by the fact that the maximal values of B; have been missed 
in this dataset, due to the lack of observations between 
JD = 2455545 and JD = 2455730. We also noticed that 
B; displays almost sinusoidal temporal variations (Fig. Q]) 
that would make the large-scale field distribution roughly 
consistent with a simple dipolar configuration (Landstreet 
& Mathys 2000). 

The chromospheric activity measurements derived from 
Call K, Call IRT and Ha, as well as their uncertainties, are 
also listed in Table[T]and over-plotted with the B; estimates 
in Fig. [T] Chromospheric emission, especially for Call K, 
follows the same general variations as B/ . This good agree- 
ment provides us with another hint that the surface mag- 
netic geometry of /3 Ceti is topologically rather simple. If B; 
provides us with a selective tracer of the largest-scale mag- 
netic structures (owing to the mutual cancellation of polar- 
ized signatures of nearby magnetic regions with opposite 



polarities), the chromospheric emission is sensitive to the 
field strength only and is therefore unaffected by this spatial 
filtering effect. The good temporal agreement between B; 
and chromospheric proxies is therefore indicative of a lack 
of smaller-scale magnetic regions outside rotational phases 
of maximum B; , in striking difference with very active cool 
stars on the main sequence for which B; and chromospheric 
emission seem to be mostly uncorrelated (Morgcnthaler et 
al. 2012). 

Temporal variations of the radial velocity are shown in 
the lower-right panel of Fig. [TJ displaying values in the 
interval 13.24 - 13.47 km/s, with relatively good correla- 
tion with other activity indicators. In early measurements, 
Frost (1923) reports RV = 13.5 km/s, while Buscombc 
& Kennedy (1968) obtained a slightly higher value of 
14.1 ± 0.33 km/s. A later estimate of 13.32 ± 0.05 km/s 
was proposed by Massarotti et al. (2008), in good agree- 
ment with our series of values. 



3.2. Zeeman Doppler Imaging 

The first step towards the ZDI reconstruction of the mag- 
netic topology of /? Ceti is the determination of the rota- 
tional period of the star. We conducted the period search 
from the full dataset (June 2010 - January 2012), follow- 
ing the approach of Petit et al. (2002). Two hundred ZDI 
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models were computed (following the method detailed in 
Section 2), assuming a different rotation period of the star 
for each model, with a uniform sampling of rotation periods 
in the 100 to 300 days interval. Forcing a constant informa- 
tion content (i.e. a same average magnetic field strength) 
in all 200 magnetic topologies, we compared them to the 
observed Stokes V profiles and analyzed the \ 2 variations 
over the period span (Fig. [2]). Our best magnetic model, 
identified by the lowest \ 2 value, suggests a possible period 
of 215 days, which we use for the final magnetic map. A 
slightly shorter period of 199 days was derived by Jordan & 
Montesinos (1991) using the value of wsini = 3 km/s given 
by Gray (1982). We applied the same maximum-entropy 
method to look for surface differential rotation (Petit et al. 
2002), but did not obtain any conclusive result. The rota- 
tional phase of every Stokes V profile (Table [I} is computed 
according to the rotation period of 215 d, with a phase ori- 
gin at JD = 2454101.5. 



5,5 



5,0- 



O 4,0- 
3,5- 




3,0 -I 1 1 1 1 1 1 1 1 1 

100 150 200 250 300 

Period (d) 

Fig. 2: Periodogram obtained using the whole dataset June 

2010 - January 2012. 

It was found from the literature that j3 Ceti has a rota- 
tional velocity in the interval v sin i = 3 — 4 km/s (3.3 ± 0.8 
km/s from Smith & Dominy (1979); 4 ± 1 km/s from Fekcl 
(1997); 3 ± 1 km/s from Carney et al. (2008)). Varying the 
value of v sin i over this range in our ZDI models did not 
significantly affect the resulting maps and \ 2 values. In the 
rest of this paper, we therefore chose an intermediate value 
vsini = 3.5 km/s. The inclination angle i of the stellar spin 
axis is the last ZDI input parameter. A value of i = 60° was 
proposed by Sanz-Forcada et al. (2002), following parame- 
ters from Gray (1989). Combining our value of usini, our 
estimate of the stellar radius (see Section 4) and the rota- 
tional period derived from ZDI, we also adopt here i = 60° 
in the tomographic models. 

In this study we present two magnetic maps for j3 Ceti, 
choosing to split our time-series into two subsets (one for 
June 2010 - December 2010 and the second one for June 

2011 - January 2012, see Fig. [3] and Fig. 01 respectively). 
This strategy is meant to highlight possible changes in the 



large-scale surface field between the two epochs (if any). 
The subsets obtained for 2010 and 2011/2012 cover a sim- 
ilar time-span and contain 25 and 13 spectra, respectively. 
According to the rotation period of 215 days, both subsets 
provide us with sufficiently dense phase sampling to recon- 
struct a magnetic map, in spite of a lack of observations in 
the second dataset at phases higher than 0.7. 

Both ZDI maps (illustrated in Fig. [3] and Fig. 0]) cor- 
respond to magnetic models with \ 2 = 2.1. The Stokes V 
data are therefore not fit at the noise level, mainly because 
our simplified magnetic model does not offer the possibil- 
ity to generate asymmetric Zeeman signatures, while the 
LSD profiles of /3 Ceti clearly display such asymmetry (with 
blue lobes generally deeper than red lobes). Such abnormal 
Stokes V profile shapes were previously observed on sev- 
eral active stars (e.g. Auriere et al. 2008) and most likely 
originate in vertical gradients of velocity and magnetic field 
(Lopez Ariste 2002). Previous attempts to get a better fit 
by implementing an ad hoc asymmetry in the model (Petit 
et al. 2005) did not significantly modify the resulting mag- 
netic topology. 

Most properties of the two magnetic maps are very alike. 
At both epochs, the magnetic geometry is dominated by 
the poloidal component of the magnetic field, which con- 
tains about 96% of the reconstructed magnetic energy. Most 
magnetic energy of the poloidal component is also stored in 
spherical harmonics modes with I = 1, with 83% and 85% 
(for the first and the second map, respectively) of the re- 
constructed poloidal energy This dominant contribution of 
low-order modes unveils a very simple topology, dominated 
by a dipolar configuration (the contribution of which is es- 
pecially visible in the radial field component). The dipole is 
almost aligned on the spin axis, leading to a majority of the 
magnetic energy in modes with m = 0. The visible pole of 
the dipole features a positive field polarity and a polar field 
strength close to 20 G (see Table [2]), in agreement with 
positive B; values recorded for (3 Ceti. We also report in 
Table [2] a limited shift in rotational phase of the magnetic 
pole between the two maps. This apparent phase drift may 
be linked to the phase gap in the 2011/2012 dataset close 
to the positive pole (at phases 0.7 and above), which may 
have induced a less accurate longitudinal positioning of the 
pole. 



Table 2: Dipole characteristics of j3 Ceti. 





rotational phase 


latitude of 


magnetic 


Epoch 


of maximal 


maximal 


dipole 




field strength 


field strength 


strength 


2010 


0.1 


52° 


20.5 G 


2011/2012 


0.0 


50° 


20.6 G 



To check again the good consistency of both magnetic 
topologies, we computed another ZDI map (not shown 
in this paper) combining all the data from June 2010 to 
January 2012. This global magnetic model is obtained with 
a x 2 of 2.6, slightly larger than the \ 2 °f the two separate 
maps. If the slightly worse Stokes V fit suggests a (lim- 
ited) temporal evolution of the field structure, the resulting 
topology is still in agreement with both topologies obtained 
from individual subsets, which indicates that the temporal 
changes have left the dipole mostly untouched and are more 
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rot. phase 



Fig. 3: /3 Ccti in the period June 2010 - December 2010. Left panel: Normalized Stokes V profiles - observed profiles 
(black lines); synthetic profiles (red lines); zero level (dashed lines). All profiles are shifted vertically for display purposes. 
The rotational phases of observations are indicated in the right part of the plot and the error bars are on the left of 
each profile. Right panel: The magnetic map of (3 Ceti. The three panels illustrate the field components in spherical 
coordinates (from top to bottom - radial, azimuthal, meridional). The magnetic field strength is expressed in gauss. The 
vertical ticks on top of the radial map show the phases of observations. 




rot. phase 



Fig. 4: Same as Fig. [3J for the period June 2011 - January 2012. 
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likely related to smaller-scale features of the magnetic ge- 
ometry. 

4. Mass, evolutionary status, and surface 
abundances of /3 Ceti 

The position of (3 Ceti in the Hertzsprung-Russell diagram 
is shown in Fig. [5] We adopt the effective temperature of 
4797 K from Massarotti et al. (2008) and the luminosity 
is computed using the stellar parallax from the New 
Reduction Hipparcos catalogue by van Leeuwen (2007), 
the V magnitude from 1997 Hipparcos catalogue, and the 
bolometric correction from Flower (1996; BC = -0.403); 
the error bar on luminosity reflects only the uncertainty 
on the parallax. The corresponding stellar radius obtained 
with the Stcfan-Boltzmann law is 18 Rq. 




7000 6000 5000 

Teff 



Fig. 5: Position of /3 Ceti in the Hertzsprung-Russel diagram 
using the effective temperature by Massarotti et al. (2008) 
and the luminosity derived from Hipparcos data (see text 
for details). Standard evolution tracks at solar metallicity 
by Charbonnel & Lagarde (2010) are shown in black for 
different initial masses as indicated by labels, with solid and 
short-dash lines indicating the phases of evolution before 
and after ignition of central helium-burning, respectively 
The blue long-dash line is for the 3.5 M© rotating model 
(see text). 

We infer a mass of 3.5 M based on the standard evo- 
lution tracks at solar metallicity from Lagarde et al. (2012) 
that are plotted in Fig. [5] According to our models, the 
main sequence progenitor of f3 Ceti was a late B-type star. 
Maggio et al. (1998) give an A-type star for a progenitor. 
Our values for the radius and mass are compatible with 
those found in the literature (15 - 17 Rq and 2.8 - 3.5 Mq; 
Jordan & Montesinos 1991; Maggio et al. 1998; Schroder et 
al. 1998; Gondoin 1999; Allende Prieto & Lambert 1999; 
Massarotti et al. 2008; Berio et al. 2011). 

As seen from its position in the HRD in Fig. (5J (3 Ceti 
could be presently cither at the base of the red giant branch 



or undergoing core-helium burning. The latest option is ac- 
tually strongly favored because of the much longer duration 
of the central helium-burning phase. In addition, the carbon 
isotopic ratio 12 C/ 13 C = 19 ± 2 (Luck & Challcncr 1995; 
Tomkin et al. 1975) and the Li abundance N(Li) = 0.01 
(Luck & Challener, 1995) derived for /3 Ceti indicate that 
this star has already fully undergone the first dredge-up. 
Indeed as can be seen in Table |3] where we give these quan- 
tities for standard and rotating 3.5 M models at the lu- 
minosity and effective temperature of f3 Ceti both for the 
first RGB ascent and helium-burning options, the observa- 
tional data are in excellent agreement with the theoretical 
post-dredge-up predictions for the surface carbon isotopic 
ratio and lithium abundance of the 3.5 Mq, Zq rotating 
model shown in Fig. [5] (dashed blue track). Clearly, first as- 
cent is excluded by the observed carbon isotopic ratio since 
at this position the first dredge-up is still ongoing while it 
is fully completed when the star undergoes central-helium 
burning (note that a 3.5 M Q star does not undergo thermo- 
halinc mixing on the RGB; see Charbonnel & Lagarde 2010 
for more details). Additionally the low lithium abundance 
for such an intermediate-mass star is a robust signature 
of rotation-induced mixing when the star was on the main 
sequence (see e.g. Palacios et al. 2003). 

Note that the 3.5 Mq rotating model was computed 
with the same input physics and assumptions as Lagarde 
et al. (2012). In particular, initial velocity on the zero age 
main sequence is chosen equal to 45% of the critical velocity, 
which corresponds to the mean value of the observed veloc- 
ity distribution for cluster and field B type stars (Huang 
et al. 2010). At the location of (3 Ceti (helium-burning), 
the theoretical surface rotation velocity of the model is 
6 km/s, in agreement with the observed usim for this star 
although it corresponds to a rotational period of 160 days, 
i.e., slightly lower than the value we derived through ZDI 
in Section 3.2. 

Table 3: Theoretical predictions for surface Li abundance 
and carbon isotopic ratio at the luminosity and effective 
temperature of (3 Ceti on the first ascent of the red gi- 
ant branch and on the central helium-burning phase in the 
3.5 Mq standard and rotating models for two initial rota- 
tion velocities (50 and 140 km.s -1 ). These numbers have 
to be compared with the observational values for (3 Ceti, 
namely N(Li)= 0.01 and 12 C/ 13 C = 19±2 (see text for ref- 
erences and details). 



N(Li) N(Li) TI C/ TJ C TI C/ TJ C~ 

RGB He-burning RGB He-burning 

standard L75 L27 67 206 

rotation (50) 1.58 1.12 61 20.6 

rotation (140) 0.53 0.007 38 18.32 



5. Discussion 

The detection of a magnetic field on (3 Ceti (with a lon- 
gitudinal strength of about 8 gauss) was first reported by 
Auriere et al. (2009). We present here our extended work 
for this star. Using spcctropolarimctric data, we derived 
two magnetic maps of (3 Ceti. We reconstructed the sur- 
face magnetic field topology of the star for two sets of 
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observations (June 2010 - December 2010 and June 2011 
- January 2012) and our model suggests that the photo- 
spheric magnetic geometry is mostly stable and dominated 
by a simple dipole. Smaller-scale features of the surface field 
topology are simply out of reach of the ZDI technique, be- 
cause of its low spatial resolution at small usini values. In 
this situation, the analysis of activity indicators is a use- 
ful complementary piece of information to partially reveal 
smaller magnetic elements. For /3 Ceti, it is clear from Fig.Q] 
that the variations of the longitudinal magnetic field B; and 
chromospheric indicators correlate rather well, giving fur- 
ther support to a simple field geometry. 

The chromospheric activity of (3 Ceti is at the level ob- 
served on giant stars with P rot « 80 d (Young et al. 1989), 
while the B/ maximal value of 8 G is typical of G and K 
giants with periods shorter than about 100 days (Auriere 
et al., in preparation). With a much longer estimated pe- 
riod of 215 d, (3 Ceti seems to exhibit an abnormally high 
activity level, at least if we bear in mind the picture of a 
global dynamo action increasing with rotation. A good in- 
dicator of the efficiency of the dynamo mechanism in late- 
type stars is the Rossby number Ro, which is defined as the 
ratio between the rotational period P ro t and the convectivc 
turnover time r c : Ro = P ro t/T c (Parker 1979; Durney & 
Latour 1978). For our 3.5 M model at the helium-burning 
position of the HRD, the convective turnover time at Hp/2 
above the base of the convective envelope is 171 days. Using 
our value of the rotation period, we obtain Ro = 1.26. It 
is generally accepted for cool main sequence stars that a 
sustained a-ui dynamo operation is observed for stars with 
Rossby numbers of the order of Ro < 1. In this respect 
again, (3 Ceti deviates from the expected trend since its 
magnetic level seems unusually high compared to its slow 
rotation. 

Since the magnetic behavior of (3 Ceti seems to be partly 
at odds with the global dynamo framework, we propose to 
explore the alternate option of a fossil magnetic field that 
may be inherited from the stable field of an Ap star, con- 
sidering that the mass estimated for f3 Ceti (3.5 M©) corre- 
sponds to a late-B/early-A star on the main sequence. This 
hypothesis was already proposed for EK Eridani (Auriere 
et al. 2008; 2011), another abnormally active cool giant 
sharing with j3 Ceti a magnetic topology almost purely 
poloidal and dominated by a dipole (in the relevant inclina- 
tion range). Assuming that this fossil field is buried in the 
radiative zone, one may question if it could permeate the 
deep convective zone of a red giant, up to the photosphere. 
Such situation has been simulated in the case of a 1 Mq on 
the main sequence, for which the convective envelope occu- 
pies about 30% of the radius. Strugarek et al. (2011; 2012) 
show that for both an aligned or inclined magnetic dipole, 
the magnetic field is able to expand from the radiative in- 
terior to the convective envelope. The interaction between 
the pre-existing magnetic field and thermal convection has 
also been studied with numerical simulations in the solar 
context. For example, the transition between regimes when 
the imposed magnetic field is strong (magnetoconvection) 
and when the magnetic field is weak (dynamo) is investi- 
gated by Cattaneo et al. (2003). 

The simple field topology of /3 Ceti is reminiscent of 
the mainly dipolar field configurations observed on main- 
sequence Ap/Bp stars (see e.g. Lueftinger et al. 2010). 
Following further the fossil field hypothesis, we can derive 
a rough estimate of the surface field strength of the main 



sequence progenitor of /3 Ceti. From our 3.5 M© model, we 
get that the radius of (3 Ceti must have been R(zams) = 
2.01 Rq when the star arrived on the zero age main se- 
quence. Assuming that the field measured for f3 Ceti simply 
differs from the main-sequence field by the dilution effect of 
its larger radius, we can use the simple formula of Stepien 
(1993) for estimating the magnetic dipole strength of its 
main sequence parent star - B(MS) = B[R/R(MS)] 2 , 
where B and R are the present dipole strength and radius of 
the star and B(MS) and R(MS) are the dipole strength and 
radius on the main sequence. In the hypothesis of conserva- 
tion of magnetic flux, we derive a magnetic dipole strength 
of B(AfS) 1650 G on the main sequence, which is con- 
sistent with typical field strength measured on Ap stars 
(Landstreet & Mathys 2000; Auriere et al. 2007). 

In the assumption of /3 Ceti being an Ap star descen- 
dant, we have finally computed an evolutionary model of 
slow rotator with magnetic braking after the turnoff. For 
Vina = 50 km/s, we did not find any difference for Li abun- 
dance or 12 C/ 13 C using models with and without magnetic 
braking after the turnoff. The numbers presented in Table |3] 
are therefore still relevant under this assumption, leaving 
unchanged our conclusions about the evolutionary state of 
P Ceti. 



6. Conclusions 

We have collected spectropolarimetric data for the cool gi- 
ant (3 Ceti in the period June 2010 - January 2012. The 
temporal modulation of polarized signatures is consistent 
with a rotation period of 215 d. Using the spectropolarimet- 
ric dataset, we were able to reconstruct two magnetic maps 
(taken one year apart) of the surface magnetic field topol- 
ogy of the star employing the Zeeman Dopplcr Imaging 
technique. Both magnetic maps display very similar mag- 
netic geometries, with a global magnetic field dominated by 
a dipolar configuration, with a dipole axis almost aligned 
on the spin axis. The polar field strength is close to 20 G 
for both maps. 

The longitudinal magnetic field Bf has significant vari- 
ations in the interval 0.1 - 8.2 G and remains of positive 
polarity for the whole observational period. Sinusoidal vari- 
ations of B; are consistent with a global dipole field con- 
figuration. The behavior of the line activity indicators Ha, 
Call K and Call IR and the radial velocity correlate rather 
well with the B; variations, giving further support to the 
dipolar model produced by ZDI. 

Based on a comparison with recent stellar models (com- 
puted with the same assumptions and input physics as in 
Charbonnel & Lagardc 2010), we derive a mass of 3.5 M Q 
for (3 Ceti and propose that this star is already in the 
central-helium burning phase rather than starting the first 
ascent of the red giant branch; the evolutionary status is 
also supported by comparison between observed and theo- 
retical values for the lithium abundance and carbon isotopic 
ratio at the stellar surface. Based on our evolutionary mod- 
els, we derive a convective turnover time equal to t c = 171 
days and we estimate the Rossby number of (3 Ceti to be 
Ro = 1.26, suggesting that a global dynamo alone may not 
be able to account for the high activity of the star. 

We derived here four useful clues to better understand 
the magnetic nature of f3 Ceti with (a) a mass of 3.5 M Q , 
consistent with typical masses of chemically-peculiar stars; 
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(b) a long rotation period of 215 d that seems unable to ac- 
count for the high activity recorded for this object through 
dynamo action alone; (c) a simple and stable magnetic 
topology consistent with dominantly dipolar field geome- 
tries observed in Ap/Bp stars; (d) a polar field strength of 
about 20 G which, when back-extrapolated to the main se- 
quence, is consistent with typical field strengths of Ap/Bp 
stars. Considering all this information together, we propose 
that the magnetism of j3 Ceti may be (at least partly) of 
fossil origin and inherited from a main sequence Ap/Bp 
star. 

The present study, together with the studies of EK Eri 
(Auriere et al. 2008, 2011) and 14 Ceti (Auriere et al. 2012) 
may give a first view on the magnetic field evolution of 
an Ap star at different stages after the main sequence, 
Hertzsprung gap, base of the red giant branch and the he- 
lium burning phase. 
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